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Abstract 

In this introductory review article, we explore the special relativistic equations 
of particle motions and the consequent derivation of Einstein’s famous formula 
E = me?. Next, we study the special relativistic electromagnetic field equations 
and generalizations of Lorentz equations of motion for charged particles. We then 
introduce the special relativistic gravitational field as a symmetric second order 
tensor field. Particle motions in the presence of static gravity are explored which 
could be used to study planetary dynamics, revealing perihelion shifts. Next, we 
investigate the system of consisting of pressureless plasmas and neutral perfect 
fluids coupled to the gravitational field. In that arena, we derive the relativistic 
Euler equation. Einally, we investigate the relativistic dynamics of a perfect fluid 
plasma and extensions to viscous flow and derive the relativistic Navier-Stokes 
equation. 
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1 Introduction 


A century has elapsed since the momentous discovery of the special theory of relativity 
by Einstein pp. This theory has inspired subsequent pursuit of the general theory of 
relativity as a novel theory of gravitation P]. It has also led to the understanding of 
the relativistic wave equation of an electron PJ. In the last five decades, relativity 
has also inspired the gauge field theories of subatomic particle interactions which have 
spectacular experimental confirmations However, applications of relativity theory 
into applied mathematical problems are almost non-existent. We venture to write this 
review article mainly to attract the academic attentions of applied mathematicians to 
this fascinating branch of modern theoretical science. 

There exist well established special relativistic particle mechanics, relativistic 
fluid mechanics and electrodynamics. In the special theory of relativity, spacetime 
is assumed to be a flat differentiable manifold. Since Einstein’s gravitational theory 
involves a curved pseudo-Riemannian manifold, special relativistic dynamics of various 
macroscopic systems usually must omit gravitation. This is generally valid as the 
gravitational field is extremely weak in comparison to the other force fields involved. 
However, a linearized version of Einstein’s theory of gravitation can be incorporated 
within the framework of special relativity. Our present review article aims at such a 
treatment of various dynamical systems. It is hoped that such a review will serve as 
a useful introduction to the field for practitioners of non relativistic fluid mechanics 
as well as those wishing to study weak-field gravitating particles, fluids and plasmas 
such as are found in various astrophysical systems. As mentioned above, it is written 
with applied mathematicians as the primary intended audience. However, we hope 
this review will be useful to the wider audience as well. No previous knowledge of 
special relativity or Minkowski tensors is assumed. 

The electrodynamics and mechanics of relativistic continuous media, coupled 
to gravitation, has many interesting applications in astrophysics (for examples, see 
inj) and other areas of general relativity and high-energy physics [0]. As well, the 
intrested reader is refered to the books [7j and references therein. 

In section-!!, notations for vectors and tensors in three and four dimensions 
are laid out in a leisurely fashion. In the third section, particle mechanics (mainly 
special relativistic) is discussed in a nut-shell. We do derive Einstein’s famous equation 
E = mcfl in this section. 

In section-IV, we discuss Maxwell’s equations of electromagnetic fields. These 
equations are known to be already relativistic! However, Lorentz’s equation of motion 
requires some minor modification for the relativistic conversion. 

In section-V, a special relativistic version of the gravitational field equations 
is investigated. It involves a symmetric second order tensor field in spacetime to 
represent the gravitational force. We then couple the gravitational field with: (i) 
an incoherent dust, (ii) an electrically charged dust (pressureless plasma), and (iii) a 
perfect fluid. In the following section we specialize furthermore to static gravitational 
fields. Planetary motions are investigated in the static field and an approximation to 
the famous perihelion shift is derived. Moreover, we also explore a perfect fluid in the 
presence of static external gravitation, deriving the relativistic Euler equation in the 
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process. 


In the last section, we generalize preceding investigations to more complicated 
materials (plasmas with pressure and viscosity) and derive the relativistic Navier- 
Stokes equation. Furthermore, generalizations to curvilinear coordinates and or¬ 
thonormal (or physical) components are also touched upon. 

In deformable media, a stress tensor, aij{x^, x^) = aji{x^, x‘^, x^) indicate 




dxi 


equilibrium conditions by the satisfaction of the equations 
the special relativistic generalization, an energy-momentum-stress tensor T^,y(x^, x^, x^) 

Tpfj,{x^, x^, x^, x'^) is introduced. Moreover, dynamical equilibrium of defovmahlehod- 


= 0. In 
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ies or fluids are characterized by conditions Ylt=i 


(x^,x‘^,x^,x‘^) p, rpl 

^—- = U. these equa¬ 
tions are of utmost importnace throughout the paper in deriving relativistic continuity 
equations and relativistic equations for streamlines. 


2 Notations and conventions 


The three-dimensional physical space is assumed to be Euclidean. A typical point in 
this space is denoted by x := (x^, x^, x^) e We use mostly a Cartesian coordinate 
system. 

The spacetime continuum, M 4 , is assumed to be a flat differentiable mani¬ 
fold admitting Minkowskian coordinate systems [H]|l()|. (These are generalizations of 
Cartesian coordinates.) Relative to a Minkowski coordinate system, an idealized point 
event in M 4 can be mapped uniquely into the point x := (x^, x^, x^, x^) G M^. Here, 
(x^, x^, x^) indicates the spatial coordinates in the Euclidean whereas x^ is the 
speed of light times the time coordinate, i.e. x^ = ct. 

Roman indices are used for three-dimensional spatial components of vectors 
and tensors (the components in the Euclidean M^). Greek indices are used for the 
four dimensional spacetime components of vectors and tensors. Einstein’s summation 
convention is followed for both Roman and Greek indices. 

We shall furnish some simple examples, three-dimensional vectors are denoted 
by a bold letter. For example, the vector v has components u*; i G {1, 2, 3}. The 
Kronecker delta: <5^-, 5ij, 6'^^ are all entries of the 3x3 unit matrix [/] 3 x 3 - Therefore, 



I for i = j, 
0 for i 7 ^ j. 


( 1 ) 


The inner product of two (three-dimensional) vectors v, w and the length | v| are given 
by 


3 3 

V • w =6ijV^w^ =: dijV^w^ 

i=l j=l 


v| := -I- dijV^v^ = a/(u^)2 -|- (u^)^ -|- (u^)^. 


( 2 a) 

( 2 b) 
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The totally anti-symmetric numerical (oriented) tensor is defined by the compo¬ 
nents: 

{ -|-1 for {ijk) an even permutation of (123), 

— 1 for {ijk) an odd permutation of (123), (3) 

0 otherwise. 

We can express some familiar vector calculus notions with the help of eijk- For 
example. 


eijkv^w^ = (v X wr)*, (4a) 

dBUx) . 

Q^k = - (V X B) . (4b) 

In the four-dimensional Minkowski spacetime, the metric tensor components are fur¬ 
nished by 

{ 6ij for /i = z and v = j, 

-1 for/i = z/= 4, (5) 

0 otherwise. 


[d^’'] = [d^^u] ^, (6a) 

d^’^d.x =5^1 = (6b) 

For the four-dimensional vector components, the lowering of indices is accomplished 
by: 

^ (T) 

so that it follows: 

Ui — dii/U — SijVp ^ 

U4 =d^i,u'^ = —u^j ( 8 ) 

The four-dimensional inner product between two vectors is provided by: 


=6ija^V — a^b'^, (9a) 

d^ua^a^ =(o^)^ + (a^)^ + (a^)^ - (a^)^- (9b) 

Vectors are characterized as timelike, spacelike and null by: 

dafiv’^yd <0 for a timelike vector, (10a) 

dajsv^^yd >0 for a spacelike vector, (10b) 

dapv'^vd =0 for a null vector. (10c) 
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3 Particle mechanics (Newtonian and Relativistic) 

We shall start with a very brief review of Newtonian mechanics. For the sake of 
simplicity, we restrict ourselves to the case of a single point particle with mass m > 0 . 
Let the parameterized motion curve by given by 

x^ = X\t) ( 11 ) 

where x* are Cartesian coordinates of the Euclidean space E 3 and t is the time variable. 
Let the three components of the force vector be given by P{t, x, v), which are functions 
of seven real variables. The components u* represent the velocity variables. 

Newton’s equations of motion for a single particle are provided by the well 
known equations: 


m 


dfl 




These equations imply that 


d , dXHt)dX^ 

— -mdii ;-;— 

dt 2 dt dt 




( 12 ) 


(13) 


The above equations are physically interpreted as “the rate of increase of kinetic 
energy of the particle is equal to the rate of work performed by the external force”. 

We note that Newton’s equations of motion m remain unchanged in form 
(or covariant) under the coordinate transformations: 

x“ =c“ + rlx’', (14a) 

[Rhx3--=[r%], (14b) 

[RV [42] = [4]3x 3 ■ (14c) 

The equation (Ei defines an orthogonal matrix The set of transformations 

(fn^ constitute the six parameter group XC>(3;M), the isometry group of Euclidean 
three-space E 3 . 

Consider another transformation, namely a special Galilean transformation: 


x^ =x^ — v^t, 
- 22-3 3 

- Jb Jb - bU ^ 

i =t. 


(15) 


The new (hatted) frame is moving with constant velocity along the x^-axis relative 
to the old frame. The Newtonian motion laws m remain unchanged in form by the 
Galilean transformations m- 

In the nineteenth century, Michelson and Morley performed some sophisti¬ 
cated experiments regarding light propagation in vacuum [2j. The startling outcome 
of their results was that the speed of light does not change due to any (constant) 
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motion of either the source or observer. Newton’s ideas of absolute space and abso¬ 
lute time (inherent in ( 1121 ) and m) are incompatible with Michelson and Morley’s 
experimental findings jHI- In 1905 Einstein solved this puzzle by the revolutionary 
ideas that space and time are relative in regards to any motion However, a com¬ 
bined spacetime continuum, M 4 , is still absolute. The appropriate generalization of 
the three-dimensional Cartesian coordinates are the four-dimensional Minkowskian 
coordinates i m- Moreover, the correct generalization of the transformations in 


(fn^ are furnished by (figure 1 ): 

=c" + (16a) 

[L]4x4 := [l^p] , (16b) 

[T)] 4 x 4 = [V] , (16c) 

[Lf [D] [L] = [D] , (16d) 

KH := [D]-^ = [D] = [V] . (16e) 



Figure 1: The inhomogeneous Lorentz transformation in spacetime. 

A typical example of the above transformation (with = ct) is provided by: 
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We make the following comments on the above transformations: 

(i) The equation (HZD is called the Lorentz (or “boost”) transformation m [im. 

(ii) It is the correct generalization of the Galilean transformation m for a moving 
frame. 

(iii) The above transformation implies contraction of length measurements in the 
moving frame. 

(iv) It can bring about the “slowing” of time measurements (time dilation) for a 
moving observer. 

(v) The speed, |u^| of the moving observer must be strictly less than c, the speed of 
light. 

The set of transformations in (llhal) constitutes a continuous group known as 
the inhomogeneous Lorentz group or the Poincare group. It is a ten parameter group 
denoted by XC>(3,1; M). 

It follows from (CMI) that det [L] = ±1. Therefore, the inverse matrix exists 
and is denoted by 


W] 4X4 W- 

= 6“^. (18) 

Now we shall define Minkowskian tensor fields in the flat spacetime manifold M 4 . 
These are defined by the transformation properties Q, lEI, H]: 

Here, the coefficients are defined by (jlfia|) and m- The tensor fields in (dni) 

are of order r + s, where r is the contravariant order and s is the covariant order. 
Note that the right hand side of (HU condenses a sum of 4”+^ terms! Moreover, these 
tensor fields are assumed to be twice continuously differentiable. The restriction of 
the tensor fields on a parameterized curve in M 4 satisfies the transformation rules: 




( 20 ) 


As simple examples, we consider the numerical second order tensors in (jlbcl) 

and JHil) respectively. By the rules m we deduce that 


doij^ — 

d^d =/« (21) 

These special tensor components retain their numerical values under the rules (PH) . 

Now, let us consider a 0 + 0 order or scalar field W{x) which is twice differ¬ 
entiable. Further, let it satisfy the wave equation: 




- q2 q2 q2 

_{dx^)^ {dx"^)^ {dx^)^ 
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1 52 ■ 

c2 {dtf_ 


W{..) = 0. 


( 22 ) 


UW{x) := 









We can prove from (which in this case reads W{x) = W{x)), (pil) . (jlhal) and the 
chain rule of differentiation that (El implies: 


UW{x) = 


d^W(x) 

dx°^dxP 


q2 q2 q2 

{dx^)^ {dx‘^)‘^ {dx^)‘^ 


1 92 


W{..) = 0. 


(23) 


Since the speeds of wave propagation in El and El are both c, the speed of light, 
we conclude that the speed of such a wave remains invariant under any motion of an 
observer characterized by (El- 

Now we discuss another important feature of the tensor field in El (or (El)- 
The four-dimensional tensor field equation 


r^"-]ri...M.(^) = o (24) 

hold if and only if the transformed components satisfy 

= (25) 

This statement physically signifies that a natural law expressible by the vanishing of 
a tensor field remains unaltered in any rotated, reflected or moving frame. The main 
mathematical postulate of the special theory of relativity is that the natural laws must 
be expressed as tensor field equations in spacetime. 

Now we shall study the relativistic particle mechanics. For that purpose we 
have to introduce the exact definition of a parameterized curve in M 4 and its physical 
interpretation. It is easier to consider the corresponding parameterized curve in the 
coordinate space M^. See hgure 2 . 



Figure 2: A parametrized curve in spacetime. 


Let a differentiable parameterized curve into the Minkowskian coordinate 
space be characterized be characterized by 


=T’#^(n), (26) 

til <ti < U2 ■ 
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Here, the functions ^ are assumed to be continuously twice differentiable. Suppose 
that the curve represents physically the history of an idealized “point” clock moving 
in spacetime. The proper (or actual) time flow of the clock along the motion curve is 
given by i, uni 


S = s*{u) := - 

c 


dX"H"°‘{w) dX# P{w) 
~dai3 -;-:- dw. 


dw 


dw 


Here, we have tacitly assumed that 


dX#f^{u) dS"^(u) 
0 < A / daf} ; ; = C- 


du 


du 


du 


(27) 


(28) 


for all n G [tti, U 2 ]- The suppositions above are based on physical principles. Consider 
for example a free particle in the spacetime. The history curve, or world line, of such 
a particle would be a straight line with some slope. An observer “traveling” with this 
particle, feeling no acceleration can justifiably claim that he or she is not moving but 
it is the surroundings which are moving. Therefore, to the observer, he or she is not 
moving in space but only in time. Generalizing this, a straight time-like world line 
therefore represents the appropriate time axis for an observer on this trajectory. 

Physically, JSHl) implies that the curve is timelike and the actual speed along 
the curve, as measured by an observer whose time axis is given by the axis in figure 
2 , is always less than the speed of light. 

The integral (EH), which defines the proper time, s, is invariant (or scalar) 
with respect to transformations (llbal) . Moreover, the integral is invariant under 
any smooth reparameterization of the curve characterized by 


=X*>^{u)=Xf^{y), 


y =yiu), 


dY{u) 

du 


/ 0 , 


Y{ui) :=yi <y <y 2 -= Y{u 2 ). 


(29) 


Chosing the parameter y = s, the proper time, we obtain from dUI), (UHl) and (Enj 
that 


Here, 

curve. 


dX°^{s) 

ds 


xt^ =X*'^{u) =: T^(s), Si = 0 < s < S 2 , 

1 f 


s =S{s) = 


'0 


dX^iw) dXh{w) 

—dotfj -;-;- dw, 


dw 


dw 


1 = 


= -A -d, 


^afd' 


dX^{s) dXP{s) 


dS{s) 
ds 

, dX^(s)dXl^{s) 2 
dap -= -c^ < 0 . 


ds 


ds 


(30a) 

(30b) 

(30c) 

(30d) 


ds ds 

are the four components of the relativistic velocity along the motion 


10 















In case the motion curve = X^{s) is continuously twice differentiable, 
which we shall always assume, the differentiation of (I.SOdl) yields 


, dX^{s) d?X^{s) _ 

- = "■ 

Therefore, in the proper time parametrization, the four-acceleration is always (Minkowskian) 
orthogonal to the four-velocity! 

In case we reparameterize the curve by choosing y = t, the usual or coordinate 
time, we derive from ((77|). (UHl) and that 



xU — =:X^{t), X* = X^{t), = T {t) := ct, 

tti 




dX'^iw) dX^(w) 

—dag —1 - 1 - dw, 

dw dw 


dS{t) 1 


dt 


= -\ -d. 


dX^^it) dX^{t) 


0/3" 


dt dt 


> 0 , 


'dS{t)' 

1 

’ dX\t)dX\t) 

dt 

c2 

dt dt 


Recall that the Newtonian 3-velocity variables in m are 

dX\t) 


=V\t) -.= 


dt ’ 


|v|2 =\Y{t)f = 6,jV\t)V^{t)>0. 
Thus, we deduce from (1^ and dsi that 

|V(f)|2 


^ds{ty^ 


dt 


=1 - 


r .2 ’ 


0 < 


dS{t) 

dt 


= + \ 1 - 


|V(t)|^ 


r.2 - 


< 1 . 


(32a) 

(32b) 

(32c) 

(32d) 


(33) 


(34) 


Now we investigate the relationships among Newtonian 3-velocity components 
with the correcponding relativistic 4-velocity components. Using the chain rule of 
differentiation, we obtain from (I3()al I3nbl I3()cll . (|Sli|32blE13), and dni) that 


d'T{t) dS{t)dXys) 


dt 

dt 

ds ’ 

dXys) 

vHt) 

ds 


|vd)P' 



dX^{s) 

1 

dX‘^{f) c 

ds 

'ds{ty 

dt 

dt 1^ 


(35a) 

(35b) 

(35c) 
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Now we shall generalize the Newtonian equations of motion m into the 
relativistic arena. We postulate, as a generalization of m, the special relativistic 
equations of motion as 0, cni: 


m 


(PX^{s) 

ds"^ 




(36) 


Here, (x) := (x^, x^, x^, x"^) and (u) := (u^, u^, u^, u^). Each of the four components 
is a function of eight real variables and represents physically the relativistic force. 
(We have tacitly assumed that the mass , m > 0, remains unchanged along the motion 
curve in spacetime.) It follows from (1311) . (I35bl I35cl and (IHHl) that 






ds 
dX^{s) 
ds 


ds 


-1 




vKt) 


(37a) 

(37b) 


Now we shall compare the relativistic equations (EH) with the Newtonian equations 
m and m- Let us digress slightly. Suppose that we have a differentiable function 
/ defined along the motion curve. By the equations (I3()all . (I32all and EH and the 
chain rule, we get 


m :=/ {S{t)] = fis), 
dfjs) _ 1 dj{t) 

ds _ |v(t)|2 dt 

Substituting into the first three equations in EH we deduce that 


(38a) 

(38b) 


A 

dt 


mV^{t) 



m 


d'^'^ (t) 
dt^ 





=P{.\xo 



(39a) 

(39b) 


Comparing (I39al) with (HH we conclude that the Newtonian momentum components, 

mV^{t) have to be modified into mV^{t)/iJl — in relativity. Moreover, the 

Newtonian force components relate to relativistic force components by the equations: 


r(Lx,v) 


- A / 1 


|V(t)|^ 


■F^x.u) ^ 4 . i i 

\ [t),x^=ct,u^=v^ 


(40) 
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The fourth equation in (ESI) yields, with ESSl), (1^ . EHEl) and ESI) 







(41a) 

(41b) 


Expanding (I41b|) for |V(t)|/c < 1, we obtain that 


A 

dt 


mc^ + ^m|E(t)|^ + 0 
1 


= 5., [v^P{..)], 


(42) 


Comparing the above equation with the corresponding Newtonian equation m , we 
conclude that the instantaneous energy E(|V(t)|) of the particle must be furnished by 


^(|V(t)|) = 


me 


1 - 


|VW|- 


= mP + ^|V(t)|"' + 0(^1. 


In the limit |V| ^ 0, we derive that 


E{0) = mc^ . 


(43) 


(44) 


The above equation, which is the most famous formula of modern science, reveals the 
enormous rest energy of a massive particle. (The rest energy associated with a single 
1 kilogram object is approximately 9 x 10^® Joules, enough to meet New York City’s 
energy requirements for more than seven months!) 


4 Electromagnetic fields 

Maxwell’s equations of electromagnetic fields are, in Gaussian units, the following: 


■ (^ijk 


aH^(x, 

dx^ 

. 

■■A 

1 

II 

^ dEi{x,t) 

dt ’ 

(45a) 

dE\-yL, 

dx^ 

. 

=4Tra(x, 

t), 

(45b) 


t) 

=0, 


(45c) 




9E^(x, 

dx^ 

. 

dB\x, 

dt 

t) 

(45d) 


’ ^ijk 


Here, E*(..) and H®(..) stand for the electric and magnetic field components respec¬ 
tively. Also, The charge density and the current density components are denoted by 
ct(..) and ji{..) respectively. 
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Another popular system of units, which we shall not employ here, is the 
Systeme Internationale (SI). In these units. Maxwell’s equations read: 


dB^{x,t) _ 1 dEi{x,t) 

^ijk 


dx^ 

dE^{x,t) cr(x, t) 
eo ' 

dB\x,t) 


dt 




= 0 , 


k ■ 


dE^{x,t) dB\x,t) 


dx^ 


dt 


(46a) 

(46b) 

(46c) 

(46d) 


The constants cq and //q respectively represent the permittivity and permeability of 
free space with the relation c = 1/yTo/Io. 

The equations (gSi-llMl) imply that the charge current conservation equation: 


da{x,t) ^ 
dx^ dt 

must be satisfied. 

The energy density of the electromagnetic field is characterized by 
u{x,t) := ^6,, [E\..)E^{..) + B\..)B^{..)] . 

The momentum density of the electromagnetic field is provided by 

Si{x,t) := ■^eijkE^x,t)B’^{x,t) 

The corresponding vector field S(x, t) is also known as the Poynting vector. 
Maxwell’s electromagnetic stress tensor is furnished by 


(47) 


(48) 


(49) 


M,,(x,t) 


+ 


1 


Ei{x,t)Ej{x,t) - -5ij6''^Ek{x,t)Ei{x,t) 


Bi{x,t)Bj{x,t) + -6ij6'‘'-Bk{x,t)Bi{x,t) 


(50) 


The components of the Lorentz force on a charged particle (of net charge e) 
is given by 


fi{t,x,v) := e 

Note that this equation yields 


Ei{x, t) + eijk—B'"{x, t) 


(51) 


E fi{t,x,w) = eEi(x,t)v\ (52) 

indicating the well known result that the magnetic field makes no contribution to the 
rate of work. 
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Now, we shall obtain the special relativistic versions of the various electro¬ 
magnetic equations. Following Minkowski m, we define the four-dimensional electo- 
magnetic field tensor components as: 




0 

Bfix) 

-B2{x) 

Ei{x) ' 

-Bsix) 

0 

Bi{x) 

E2{x) 

B2{x) 

-Bi{x) 

0 

Efix) 

-Efix) 

-E2{x) 

-Efix) 

0 


(53) 


It should be mentioned that Minkowski first unified electric and magnetic fields by 
the definition (I5H1) of the four-dimensional anti-symmetric tensor F^^ix). The above 
is one of several possible definitions for the tensor We define the relativistic 

charge-current components by: 


Jfix) :=/(x,t), 

J^(x) :=C(j(x, t). (54) 


The Maxwell equations (l45aM5dl) . along with definitions (p|) and (El, neatly boil 
down to m 


^F^^^{x) 

dx'^ 

dFi^‘'{x) dF^^{x) dF^f^ix) 

dx^ dxf^ dx’^ 

dJ>^{x) 


dx^^ 


dvr 

= 

c 

= 0 
= 0 . 


Outside of charged matter, Maxwell’s equations are summarized by: 


dF^,y{x) 

dx^ 


+ 


dFt^^jx) _ 
dx’^ 

dFyx{x) ^ dFxi,{x) 
^x^^ dx'^ 


which together imply 


F\F^y{x) 


' 0 ^ d'^ 

_{dx^)^ (9x2)^ {dx^)‘^ 


1 ' 


Ff,u{x) = 0 . 


(55a) 

(55b) 

(55c) 


(56) 


(57) 


Since the above equations are tensor field equations in the four dimensional spacetime. 
Maxwell’s equations outside matter were already relativistic even before Einstein’s 
discovery of the special theory of relativity! Moreover, the equation El implies that 
electromagnetic waves propagate with the speed of light. Since the wave operator 
□ is a relativistic invariant (see equation dH), we can conclude that the speed of 
electromagnetic wave propagation remains unchanged under the boost transformation 
(|16a|) characterizing a moving observer. Thus, Michelson and Morley’s experimental 
puzzle is logically explained. 
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Now we shall unify energy density, momentum density and Maxwell’s stress 
tensor by defining a relativistic electromagnetic energy-momentum-stress tensor: 


:= 




1 

dvr 




—M®-^(x, t) cSj(x, t) 
c5j(x, t) u{x,t) 

dMj 


= >1^", (58a) 

(58b) 


Working out the covariant divergence, -gff-, using (I55al) . (I55bll and (I58al) yields 


dMcf 


1 


dTT 

1 


dF, 






p\ \ _|_ p f3 dF\ ^ 

"" " " ^ dx!^ 


(9x° 


=-F T\ _F 

=—F\J\ + 0. 

c 


A/3 


^ ^ ^ 

(9x" 


(59) 


Now, the relativistic Lorentz equation for a charged particle with mass m and 
charge e are taken to be (compare with (jSSJ) 


m- 


d ^X^{s) 

ds"^ 


=F%..)\.. := -Fl{x)\^^s) 


dX^{s) 

ds 


^ , , dX‘(8)dX\,) 

-^ = 0 . 


(60a) 

(60b) 


1 - 


|v(t)p 




^ dXUt) ^ , 

Fij{x, h cFi4(..)| 


=e 




J \x(t) 


(60c) 


The right hand side of (I60cjl yields the correct Lorentz force given in equation m- 
Now we shall briefly introduce the electromagnetic four-potential. According 
to the converse Poincare lemma d, the equations (I55b|) imply that there exists 
relativistic field components A^(x) of class such that 


Ffj,u{x) - ^^^Ay{x) ^^^A^{x), 


so that 


Ei{x,t) = 


dAi 1 dAi 
dx^ 


dt ’ 




dAu dAd 


dx^ dx^ 


(61) 

(62a) 

(62b) 


We note that the four-potential components A^{x) are not unique. We can make a 
(local) gauge transformation: 

dA{x) 


=A^{x) - , 


(63) 
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Here, A(x) is an arbitrary function of class C"^. We can perform a gauge transformation 
such that the function A{x) satisfies the partial differential equation 


□A(x) 


dA>^ 


Thus, from equation and dsi, we get 


(64) 


dA'i^ 

dx>^ 


dA^ 

dx>^ 


□A = 0. 


(65) 


The above condition on is called the Lorentz gauge condition. In this gauge. 
Maxwell’s equations (|55al - [HH^ reduce, by H65j) to 




d 

dx^ 


dA 


lu 


dx’' 


- = 47rJ^(x) 

dJ^{x) 


^x^^ 


= 0 . 


(66a) 

(66b) 


The relativistic inertial energy-momentum-stress tensor for the charged matter is 
given by: 


[I^’^ix)] 


■ p{x)U\x)W{x) p{x)W{x)U^{x) 1 
p{x)W{x)U^{x) p{x) (C/"^(x))^ 

d^,U^{x)U''{x) = -c\ 


(67a) 

(67b) 


Here, p{x) represents the invariant (or proper) mass density and U^{x) are components 
of the four-velocity field. Thus, p{x)W{x) are related to the three-momentum density 
components and p{x) [C/^(x)]^ is the energy density. 

Observing the similarity between Ad“^(x) in ()58bl) and X“^(x) in H67a|) , we 
define the total energy-momentum stress field for a charged material as 

r“^(x) := T“^(x) + A4“^(x) = T^“(x) + A4^"(x) = T^"(x), (68) 


Using (iM . ((sni), (EMiniEi) and (|M1) we derive that: 


dT'^f^ix) 

dxd 


pm 


dxd 


+ u' 


,d{pu^) ^ \p. 
c 


Xoi 


dxd 


Jx- 


We postulate that a physical conservation law to hold: 


QjpaP 

dxd 

Moreover, we consider the case of a charged dust, i.e. 


J^(x) = a(Q){x)U^{x). 


(69) 


(70) 


(71) 
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Here, cro(x) is the proper charge density (the charge density measured in the co-moving 
frame of the charge). The equation with (I67bll . yields 


d 

dx^ 


p(x)U^(x) 


pU^(x) 


dU^{x) 

dxf^ 


= 0 , 


c 


U\x)F\ 


(72a) 

(72b) 


The equations (I72al) stands for the continuity of the material flow, whereas the equa¬ 
tion (l7^ represents the Lorentz equation (jtiOal) for the charged dust. 

Following the discussion on particle mechanics, we define the three-dimensional 
velocity field components by: 


u* = H*(x, t) U^{x) = 


F = H*(x,t) =\ I- 


m{xj 

V(x,t)|2 


|V(x)|2 

—— 


U^{x), and J^{x) = ao{x)U‘^{x) = ccj(x). (73) 


The spatial components (( 22 EI) provide: 

..k d 


P 




1 - 


V 


dx^ 


1 - 


= O-Q 




(74) 


Using notations of three-dimensional vector calculus, the equation dZH) can be cast 
into the more familiar form: 


m 


V(x,t) 


1 - 


iv(x,bp 

— 72 — 


-Pv-V 


V(x,t) 


1 - 


iv(x,bp 

— 72 — 


'1 


r V 


> = O-Q 

E -1— X B 

A 

I 

c 


(75) 


5 Special relativistic gravitational fields 

5.1 Introduction 

The static Newtonian gravitational potential, W{x), satisfies: 


VV(x) :=d^^ 


d^lU(x) 

dx^dxi 


47 rGp(x) inside a material body, 
0 in vacuum; 


,r(x) = - O / /hi 

Jbody |X - X' 

/(x') >0, lU(x) < 0. 


■ d^xP 


Here, G is the Newtonian constant of gravitation. 


(76a) 

(76b) 

(76c) 
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The well known equations of motion of a particle of m > 0 in the external 
gravitational field W (x) are furnished by: 


dY{t) , 

= -mVTT(x)|;^(,), 

dV^{t) _ aiT(x) 
dt dx^ \x{t)' 


(77a) 

(77b) 


The natural relativistic generalization of the Newtonian potential is a four¬ 
dimensional scalar field. Such a held has been considered but the scalar held theory 
yields incorrect planetary orbits m- Next in order of complication, we should consider 
a four-dimensional vector held as in equation (j filial) . However, material sources of 
such a held repel each other (like similar electric charges). Thus, vector helds are 
ruled out as viable candidates for a relativistic gravitational theory. Next in order of 
complication is the four-dimensional second order tensor held. This type of held is 
a potential candidate as the relativistic generalization of the mass density (within a 
factor of c^) is the energy density. As discussed in the previous section, the energy 
density is a component of a symmetric rank two tensor, the energy-momentum-stress 
tensor. Therefore, it is natural to consider a potential which is a symmetric second 
rank tensor held with the energy-momentum-stress tensor as its source. The hrst-order 
approximation of Einstein’s general relativity theory of gravity yields a second rank 
theory similar to that presented below. The held equations, along with supplementary 
conditions are summarized as the following: 


dT^^'^ix) 

= 0 , 


SttG 


K := 


dx^ 

d(j)'^''{x) 


dx’' 




f 4, 2G /■ 

^ J body 


(x', x^ - |x-x'|) 3 , 

—^ d X , 


X — X' 


(78a) 

(78b) 

(78c) 

(78d) 


Physically, the components of the tensor, T^j^, represent the total density of energy, (c- 
times) momentum and stress of the source material. (A special example was furnished 
in the equation dSHl) .) The component T 44 (x), Ti 4 {x) and Tij{x) represent energy 
density, momentum density and stress (or rate of stress) density respectively. The 
equation (|78al78c|l are analogous to the electromagnetic equations (jnsi-inn3) and 
respectively. 

Now we shall express the equation of motion of a test particle of mass m > 0 
in a gravitational held. The equations, though complicated, are aesthetically pleasing. 
Firstly, it is convenient to introduce a related second order symmetric tensor held by 
the dehnition: 

gfiuix) ■■=d^y + 4>^v{x) - ]^d^ud°'^(l)ai3{x) = gut^ix), (79) 
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For weak gravitational fields, \(f)^i,{x)\ << 1 and the determinant of gf_iu{x) is close to 
— 1. Thus, there exists a unique inverse matrix defined by 

\g^'^{x)]--=[9A^)]-\ (80) 

{x)gu\{x) 

Treating formally g^u{x) as a “metric tensor”, we define the associated Christoffel 
symbols as 


a 

(5 7 


:= i^9^\x) 


dg'yx{x) dgxgjx) _ dg^j^jx) 
dxP dx"/ 


dx^ 


a 

7 /? 


(81) 


In our linear theory, the above are components of a third order four-dimensional tensor 
with forty independent components! 

We postulate that the special relativistic equations of motion, subject only to 
gravity, must follow geodesic paths of the “metric” gfj,^{x): 


d^X^js) 

ds^ 


+ 


a 

P 7 


dXf^{s) dX'^is) 


\X{s) 


ds 


ds 


= 0 . 


(82) 


Note that these equations of motion are independent of mass, a property in com¬ 
mon with the Newtonian theory. These semi-linear equations distantly resemble 
the Lorentz equations (l6Ua() and generalize Newton’s equations (IttEIi considerably. 
Geodesic equations (|8^ admit the exact first integral im 


, , dXf^(s)dX’^{s) 

9tiiy[x)\x{s) -^-^= constant. 


ds 


ds 


(83) 


Since, in the limit (j)^u{x) —> 0, the equations (IHHl) reduce to (ISOdjl and (1670 , we must 
choose 

, , dXt^{s)dX''(s) 2 / .X 

The generalization of dHl in the presence of gravitational as well as non- 
gravitational force components J-°{x,u) is given by 


d^X^{s) 

m —— = —m 

ds^ 


a 

P 7 


\X{s) 


dxf^{s)d^x'y{s) , 

\^'>'^)\P^(s)^dX{s)/ds: (85a) 


ds 


9ap{x)\X{s) ^^^^^'^ ^^ix,u)l„ = 0 . 


(85b) 


(Compare the above equations with (i36|) and (1^1 . 

With the help of the “metric” tensor g^y{x), we can define the familiar co- 
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variant derivative El and the following consequences: 

( 86 a) 

- { 1,3 } 

- { a fj } - { a 7 } 

"^aOp-fix) = 0 , VaQ^'^ix) = 0 , V„(5^ = 0 , ( 86 d) 

VadfB^ ^ 0 , VadP"' ^ 0 . ( 86 e) 

5.2 Gravitational field of an incoherent dust 


The gravitational equations (I78al - 175^ are all linear. Splitting the solution into its 
vacuum and inhomogeneous parts we express: 


(j^fiuix) 4^{i)^u{x) T ^(0)/j,i^(x)j 

dT^^Tri dcji^jJx) 

'd^(|>(I)^Jiv{x) = - 2KT^y{x), -= 0, -= 0, 


dx'' 


dx’' 


d(t)^Q){x) 


□</'(ow(a;) = 0 , = 0 , 

t ^ ^ , , , 2G f T^i.(x', x^-|x-x'|) 

(pf,u{x) =4>(o)fiu{x} + — - 1 _ n -« 

^ J test body I ^ ^1 


3a:'. 


(87a) 

(87b) 

(87c) 

(87d) 


The above is the most general solution of the partial differential equations ((THaJ-[TScI) . 
Here <i>(^i)^v{x) represents the particular solution due to the test body material with 
energy momentum stress tensor T^y{x) and <t>{Q)^p{x) represents the vacuum solution 
created by any external sources (see figure ©). 

Ignoring the internal field of the test body, we modify dzni) and (|511) to 


5 (o)/ 3 i. := + 4>(Q)^lu{x)-]^d^yd^'^(t)(Q)o,i3{x), ( 88 a) 

g{o)^u{x)U^{x)U^{x) = - ( 88 b) 

:= g(Q)^,uu'', = -c^. ( 88 c) 


For an incoherent dust coupled to gravity, the energy-momentum-stress tensor is taken 
to be 

r^''(x) := p{x)U^^{x)U^{x) + ©^"(x, p, (?i(o)7a), (89) 

Here p{x) is the proper mass density, U^{x) is the four velocity-field of the test particle 
and 0 ^'^(x, p, ^( 0 ) 70 -) indicates the interaction of the test body with the external 
gravitational field (from this point onward we will write this term as 0 ^'^(x)). 
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Figure 3: Gravitational forces on a test body due to an external source in spacetime. The 
external source’s gives rise to (/>(o)^jy while the test body’s T^i, yields 4>{i)fj,v 


The conservation equation (|78bl) along with (P|) yields: 


2d{pu’') 


, d{pu^) 

dx’^ 


+ pu’^ 


duf^ 


+ 


dx^ dx^ 


= —c 


dx^ 


+ pu^U/j, 


dx’- 


+ 


u A 


( 0 ) 


+ Un 


dx’^ 


Note that the middle terms in the above equation vanishes via 


A 

u A 


u’'u^ 


( 0 ) 


(90) 


du>^ 

dx’- 


+ 


V A 


( 0 ) 




(91) 


At this stage there are several possible choices to enforce the conservation law CT- 
We stipulate the condition; 


dx’-' 


= P 


A 

a (3 




(0) 


(92) 


as this choice yields the correct physical laws as dictated by the continuity equation 

HSI: 

A lp(^)U-{:c)] = 0. (93) 

The right hand side of © therefore represents the gravitational force components 
due to external sources. 
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With the help of equations dZl , we can express (El in a slightly more familiar 


form: 


m 




= 0 . 


(94) 


Substituting El and El into the conservation equation (I78b|) with (I86aj) and 
El we obtain 


u 


V 


du^ 

dx'^ 


u A 



u''y{0)uU^ = 0 . 


(95) 


For the streamlines, we need to solve the system of first-order ordinary differential 
equations: 

^ (a;)|A’(s)- (96) 

(We assume the validity of the Lifshitz conditions m on the right-hand-side of (E1-) 
Using the above eqnation in El we derive that 


I dX‘'{s)dX^{s) 

ds'^ 1 A j ds ds 


(97) 


Therefore, streamlines/oZ/ow “geodesics” of the external metric g{Q)^i,{x). The Christof- 
fel symbols mediate the external gravitational forces on a test particle. 


5.3 Gravitational field of a charged dust 

A fluid with vanishing pressure is know as a “dnst”. In this case, the combined 
electromagnetic and gravitational field eqnations are investigated m- Using (|55al - 
I55cf) . (ITTl) . I58al) ■ (I67al) and H87allH7cl . we obtain the pertinent equations to be: 


dvr 

dx- - c ’ 

(98a) 

Q Q p’^>‘ Q 

dx^ ' ^x^^ ' dx^ ’ 

(98b) 

dx^ “ ’ 

(98c) 

0(7)a/3(^) 

(98d) 

o' 

II 

II 

o' 

□ 

(98e) 

^4^{I)af3 — .— 2/^) [pllo^Up Qap “1“ -M-ap] , 

(98f) 

o' 

II 

1 

(98g) 

< _o 

dxd 

(98h) 
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( 99 ) 


Using (l98tD . (TTT]) and (1691) . the conservation equation (|98gD yields: 


0 = 


dxf^ dxl^ dxl^ c ^ 


dx!^ ' dx^ ' c ^ 

Imposing the condition (HD) on 0 “^(x), the above equation (inn) leads to 


+ pu^vfu^ + ^uf^Fg^ = 0 . 
OXP ^ c ^ 


( 100 ) 


Now, the equations (I 86 dll and (I 88 bl) provide: 


2UaV^°^u" =5(0)a7 




=VL°) [5(0)a7^“«1 = (-c') = 0- 


( 101 ) 


(vi°^ is the covariant derivative defined with respect to the metric g{p)^y, defined in 


^.) Therefore, (Cnil), (EDI) and (irnni) yield 

u^Fa 

dxP - c3 “ ^ 

Substituting (nni into (irnni) . we finally obtain 


( 102 ) 


puNfu^ + 


^u^F, 


a 


(5" + -^U'yU°‘ = 0. 


(103) 


Hence, the stream lines of a charged dust satisfying (EEI) must pursue trajectories 
governed by the equations of motion: 


p{x) 

(7(0) (x 


d^X^js) 

ds^ 

■Fg^ix) 


+ 


a 

(3 7 

1 


dX^{s) dX^{s) 


(0) 


ds 


ds 


K + 72£'(0)7M 


^^"(5) dX^^{s) 


ds 


ds 


dX^ji 

ds 


]x‘=‘=X°‘(s) 


= 0. (104) 


These are the modified relativistic Lorentz equations of motion in the presence of 
external gravitational fields. (Compare with the equations (( 22 bl) and 
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5.4 Gravitational field of a perfect fiuid 

Following the prescription in the previous sections, a special relativistic perfect fluid 
in a gravitational field is governed by 






= - 2 kT“^(x) 


Qrpafj 

dx^ 

a/l 

(I) 


34.“^ 


dxf^ 


:= — 2k 

= 0 , 

= 0 . 


p{x) + +p(x)g'“Q^(x) + 0 "^(x 


(105a) 

(105b) 

(105c) 

(105d) 

(105e) 


Here, p{x) is the pressure. 

The conservation equation (|105d|) implies from (llOScI) that 


^ d 

0 = .3 = u 


dxf^ dx^ - 

00a/3 


p +?) >■'’ 


+ fp+4) uNfu^ + 


d 

dxP 


0/3 

P9m 


+ 


dxP 


I , P 
~ [P+ -3 


a 

Pi 


u^v? 


( 0 ) 


(106) 


Following the previous section, reasonable physics demands that we stipulate 


500/3 

dxP 


= [P + 


P 


a 

P 7 


u^vl. 


(107) 


( 0 ) 


Note that the right-hand-side denotes the external gravitational forces on the effective 
mass density {p+ Now, the equation (I106|) yields 


u 


d 

dxP 


P + 




+ [P+^ 


= - 


d 


0/3 


(108) 


Contracting this equation with Ua and using 
equation, 


^ “ dxP ■ 

and (UnU), we deduce the continuity 


A. 

dxP - 


p+ 


d 


= ^u. 




0/3 


Using this equation in (CHHl), we derive the relativistic Euler equation 


(p+|)“M°’“” = - 


o , 




d 


ltd 


dxP 


(109) 


( 110 ) 
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On a typical stream line (given by (1961) 1 of the perfect fluid, the following 
equations of motion hold: 




ds^ 

1 

+ ^5(0)7/^ 


^ r a 1 dX^{s) dX^js) 
I /? 7 i ( 0 ) ds ds 

dX^{s) dXf^{s) 
ds ds 


^ 1 7/3 

^ r(o)J 




( 111 ) 


6 Static external gravitational fields 


The special relativistic gravitaitonal fields may easily be compared with Newtonian 
theory in this case. The relevant relativistic equations read: 


=0: 4>{o)iA{x) ^ 0, (112a) 

hoMx) ^0, = 0, (112b) 

= - 0(o)44(x), (112c) 

9(0)ij{x)= 1 + ^())(0)44(x) 6ij, (112d) 

9(0)i4{x) =0, 5 ( 0 ) 44 ( 2 :) = -1 + ^0(o)44(x). (112e) 

(We have made use of equations (l8<Sal) to derive the above.) 

The equations IHZEl reduce Laplace’s equation: 

V2</.(0)44(x) =0. (113) 

Comparing uni) with (|76al) , (|76b|) and (I87dl) we identify (()(o) 44 (x) with the Newtonian 
gravitational potential via 

<(>(0)44(x) = -^lT(x) = ^|VL(x)| > 0. (114) 


We also note that for a four-dimensional vector field r“(x), the equations (ITTMll . 
it™ and (imi) yield 


-9(o)ap{x)T°'{x)Td{x) = 


1 - 


1 “ 2^{0)4:Ai'^) 
2|W(x)| 


[T\xf- 


[tHx)Y- 


1 + 


1 + 2 ‘^(0)44 (x) 

2|W(x)| 


6ijT\x)TYx) 
5ijT\x)TYx). 

(115) 


We have implicitly assumed that < 1- The physical components of the vector 

field r"(x) are 


r{x) 

t\x) 



2|W(x)| 


T\x\ 



2|W(x)| 

r.2 


T\x). 


(116) 
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Notice that these components are just the components of T°‘{x) projected into the 
corresponding orthonormal coordinates (or frame) of the metric g{Q)^u- 


6.1 Test particle motions in external static gravitational fields 

In a previous section it was noted that any stream line of incoherent dust follows the 
“geodesic” equation (EZI). Along any of these time-like geodesics, the equation (I51|) 
implies that 

dX^^{s)dX-{s) _ 2 

-^ = 

The above condition modifies the equation (llOdl) . For the sake of consistency, we must 
alter the definition (EH) for the proper time along a time-like curve by: 




dX^{w) dX‘'{w) 


, 1 r I , , dX>^iw) dX-iw) ^ 

The equation (I.SOdl) is modified by (tra and (ITTTl) into 


(118a) 

(118b) 

(118c) 


1 + 


2|lT(x)| 


fi- 


dX\s)dX^{s) 
ds ds 


1 - 


2|lT(x)| 


I- L 


dX^(s) 
ds 


= -c^ < 0, (119) 


while the equation (jMl) is changed by (inni) into 
ds dS{t) [ 2|lT(x)| 


dt 


dt 


= Wl- 


1 + 


2|lT(x)| 


|v(t)p 


< 1 . 


( 120 ) 


I- 


The above equation reveals the time dilation along a moving dust particle in an 
external gravitational field. 

Similarly, the equations (I35bl) and (|35cB change over into: 


UHs) = 


U\s) = 


dX\s) 

ds 

dX^{s) 

ds 


V^t) 


1 - 


2|ty(x)i 


- 1 + 


2|W^(x)|\ |Vp 


1 - 


2|H/(x)| 


- 1 + 


2|H^(x)|\ |V|2 


I- 


I- 


(121a) 


(121b) 


Here, V^{t) = are the components of the “coordinate” velocity which in general 

/l I 

differ from the “measurable” velocity components, V (t) := J ^ 2 |wwi —dt \ ■ 
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Now, the geodesic equations (jSZl) are derivable from a variational principle 
HU. A free particle’s motion is governed by a purely kinetic Lagrangian, ^mSijV^V^ 
in non-relativistic mechanics. The relativistic analogue is \rn g{Q)aijU°^u^. In the 
presence of gravity, the gravitational coupling arises from lowering the index with the 
metric ^(o)^/? and therefore, using (I112dl) . (I112ej) and (jlldl) the Lagrangian (per unit 
mass) becomes 


L(x, u) := 


1 


1 + 


2|IT(x)| 


diju'-u^ - - 


1 - 


2|IT(x)| 


4\2 




=-5iju'u^ - IT(x) 


4\ 2 


{uy lT(x) 


2 SijU^uK 


( 122 ) 


The relativistic Euler-Lagrange equations are given by 


dL{x, u) 


d dL{x,u) 
ds du°‘ \.. 


= 0 . 


(123) 


The Newtonian Lagrangian for the corresponding Newtonian theory is 


LAr(x, v) 




IT(x), 


(124) 


giving rise to the equations of motion (1770. 

For small velocities, by the equations dnn, dmi) matches dHU) except for 
the term — (tt‘^)^/2. This term represents the large rest energy contribution which is 
manifest in the relativistic physics. 

It will be instructive to investigate the fourth equation of (IT^ . Since is a 
cyclic variable in the Lagrangian (HU, the corresponding equation of motion admits 
the first integral; 


dL{x, u) 
du'^ \.. 


1 + 


2IT(x) 


or. 


1 ^ 


U^[s) = const., 


8 




1 - 


2iy(x) 

^2 




= — = const, 
c 


(125a) 

(125b) 


For small velocities and weak gravitational fields, equation (I125bl) yields: 


8 = c^ + 


-|V|' 
2 ' ' 


+ IF(x) 


J|.. 




(126) 


It is not difficult to interpret the above equation. The constant 8 represents 
the conserved total energy per unit mass. The first term on the RHS is the large 
rest energy of the unit mass. The second and third terms are the usual kinetic and 
potential energies respectively. 
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6.1.1 An example from planetary motion 


Here we investigate planetary motions due to the spherically symmetric gravitational 
field of the sun. The usual potential function is furnished by: 


lT(x) = 


-GM 


^(a:i)2 + (x2)2 + (x3)2^ 


(127) 


with M > 0 being the solar mass. Employing spherical polar coordinates in space, 
the equations and (inii) yield the Lagrangian 


L{x,u) = ^ ( 1 + 


2GM\ r 2 , 2/ 0\2 I 2 ■ 2 at ip\2 

■ I [u ) +r [u ) + r sm 0(u^) 


c^r 


-^1- 


2GM\ 


c^r 


(128) 

The Euler-Lagrange equations (!I23, from dEHl) admit uniplanar motions 
characterized by: 

"" (129) 


e = Q{s)^-. 


The reduced Lagrangian is then 

= ) (l + ^) [(“')" + - ) (l - 

In this case, the conservation of energy equation (TTS^ reads: 


2GM\ 


.4\2 


c^r 


[u 


1 - 


2GM 

c^r 


U\,) = 

c 


(130) 


(131) 


The presence of another cyclic coordinate, ip, leads to the conservation of 
angular momentum: 


= const. 


du, 




c^r 


-J ^ |.. 


(132) 


Substituting (tT^ . (unu and (HSU) into (EUl), we obtain 


/ 2GM\ 2 - 2 

I, + 

We reparameterize the curve by the following: 

r =7Z(s} = 7Z((p}, 

d7Z(s) „ diZiip) h diziip) 

_ ^^7/' _ _ 

ds I" dip r2 (l + 2^^ ^ dip 

It is useful at this point to make the following coordinate transformation: 

y=-, y = Y{(j)), 

dYjip) ^ 

dip 


(133) 


(134) 


(135) 


29 




























With (IliMj) and H135I1 . the equation (ESI reduces to 


1 - 


2GMy 

—pi— 


^ _|_ 2GMy 




2GMy\ _ 

)~Wi 


(136) 


The above first-order equation can be solved by quadrature. However, to extract 
physically important effects, we differentiate the equation (EMI) to get 


y” + y = {l + 


AGMy\ 


GM 2GM 


/l2 




{{y'f + y^) 


+oi4 


GM , fGM 

''IF 


JGMY 2GM 

y + — (137) 


The equation ESIl) may be solved by the perturbative expansion: 

y — yoH yH T + --- 


(138) 


Using this expansion in ESHli we obtain 


yi GM ^ 4 (GM )3 

y = yo + -^ + ■■■ = [1 + ecos(<y 9 - a;)J H- „2 7.4 eipsm{(p - uj) + ... 




(139) 


Here, the constants of integration e and co represent the eccentricity and perihelion 
angle of the orbit respectively. Combining the first two terms in (EHi), we conclude 
that (7M 

2/0 + ^ = —[pT [1 + e cos(</i - nJ - (5uJ)] -|- (higher order), (140) 


where 


6uj := arctan 


/i2 


Aif 


(“)■ - 


GM 

~h~ 


(f -|- (higher order). 


Thus, the elliptic orbit processes and the perihelion angle changes (see figure EJ. This 
gravitational perihelion shift per revolution is given by 


Auj = Stt 


GM 


(141) 


For the planet Mercury, the above amount yields a little over 57" per century! How¬ 
ever, the full non-linear theory of general relativity predicts the observed amount of 
almost exactly 43^^ per century ^ j2j. 

^The actual observed perihelion precession of Mercury is approximately 5600” per century. When 
effects such the attraction due to the Newtonian gravitational field of the other planets are taken 
into account, along with the fact that the Earth is not an inertial frame of reference, a residual 43” 
per century persists. The origin of this residual precession was a mystery until general relativity was 
formulated in 1915 |5|. 
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Figure 4: The perihelion precession of the planets in a gravitational field. The diagram 
displays two orbits, perhaps millenia apart, of a planet around a star. The perihelion point of 
the planet has shifted by an amount AoJ, an effect not predicted by Newtonian gravity theory. 


6.2 Perfect fluid in static gravity 

In this case, the equations (Tn^ . (tTTO and (ESI) yield, for the metric tensor, 

2VF(x) 


9{0)ij 

21F(x)' 

ij 

■ 21T(x)- 

c2 

%)- 


— 0 ) 9(0)44(3^) — “ 


1 + 


1 -1 


= 9 i 0 )ix) = 


1 + 


2W(x) 


1 -1 


The corresponding non-zero Christoffel symbols are provided by: 


j k 
4 

i 4 


(0) 


(0) 


1 - 


21F(x) 


-1 


dW{x) ^ dW{x) ^ dW{x) 

—^-7- Okl „ ■ - Oij- 


4 

4 i 


(0) 


1 dW{x) 
(? dx* 


dx^ 

i 

44 


dx^ 


(0) 


1 - 


21T(x) 


-1 




(142) 

(143) 

dlT(x) 

dx^ 


The continuity equation (tTUUl) for a perfect fluid yields m 


A 

L 


d 

dx^ 


, P 1 


+ 


d 

dx^ . 


, P 1 




Using ESI), the equation (imii leads to 


■ d 

( ah\ 

d 

( «4\ 

dx^ 

(ra(o)J 

(9x^ 





(14 


p+-y u 


+ 


d 

dx'^ - 


P + ^]u 


1 - 


2VF(x) 




p 1 - 


21T(x) 


4) 


-1 


(145) 
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Substituting (I121al) and (I121bl) into (tm . we derive that 


d 

ip+ 

yb 

d 

-U 

1- 

+ 

1 _ 

dx^ 

C2/ 

^l + 2^-(l-2^)^ 

dt 



1 - 


2W{x) 


d 


^1 + 2^-( 1 - 2 ^) Hi 


dx^ 


pil- 


2W{x) 


. (146) 


Expanding the above in powers of c we can express 


d 

dx^ 




VE(x) |v| 


+ 


2c2 


d 




W{x) 


+ 


V 


2\ n 


2c2 


1 


dp 


= -r 


1 


dx^ 

In more familiar notation the above equation reads 


(147) 


^ dp 1 


dp d 

~'di ^ di 


—pV • V + V 

U.|2' 


p nE(x) - 


p IE(x) - 


1 


+ 0(74 


(148) 


The relativistic correction terms are all collected on the right hand side of (ITIHI) . 

Now we shall investigate the equations of stream lines of a perfect fluid in the 
external gravitational fluid. Equations mu) provide 


, P 

P+ -2 
& 


d?X\s) 

ds"^ 


+ 


i 

j k 


dX^{s)dX^{s) 


( 0 ) 


ds 


ds 


+ 


i 

44 


dX^{s) 


^ U(0)7i + 5(0)74 


dX'^{s)\ dX 


' i 1 dX^ dX^ 


d 

dx^ 


ds 

,bk 


ds 

1 

~ 'r- 

I- ^ 


d 

dx^ 


(0) 

67 
^^5(0) 


ux^jsy ^ 

y ds 


I- 


+ 


d 


,47 


dX^dX^ d 

,2 5 ( 0)44 


dx^ v^^(o) 

,44 


I- 


(149) 
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Using (11421) and (I14:-{ll along with (11491) . we deduce that 


, P 

P + — 


(fX^s) 1 / 21U(x) 


ds‘^ 


+ -:t 1 — 


-1 


5jkS 


iidW dW ^^dW\dX^{s)dX^{s) 


dx^ 


-S\ 


—^ _I-^-1 

dx^ ^ dx^) ds ds 


+ 


1 ^ _ 2Tu(x) y^ 

V ) dx^ \ ds J 

1 / 21U(x) \ dX^ dX^ 




dk 


^ ds ds 


J I- 
d 
(9x“ 


p 1 - 


21U(x) 


-1 


ea/c 


2TU(x) \ dX^dX^ I d 
ds ds I dx^ 


, . 21U(x) 

P(l + —^ 


-1 


(150) 


We parameterize the curve by the rule: 


A. 

ds 


dSjt) 

dt 


-1 


A 

dt 


2^ 


-(1-^ 


2W{yi)\ |v|2 


dt 


(151) 


Then, from (Csni) and (11511) the following is obtained: 


{p + 


p\(dV%t) V\t) d 


+ ^ 1 — 


dt 
2W(x) 


2 dt 


In 1 + 


- 1 - 


-1 


Sjk5' 


2W(x) 


dW i dW i dW 

dx^ ^ dx^ ^ dx^ 


21U(x)\ |v 


V^V'^ 


2H;wVU.|f 


I- 


1 + 


2TU(x) 


- 1 - 


2TU(x)\ |v 




5\ + -{l- 


2W{x) 


SkiV^V^ 


X d' 


ak 


A 

dx° 


pll- 


2W{x] 


“ 3 ' ^ + 


2W(x) 


V I 

1 dt 


, . 2W(x) 

P(l + —^ 


-1 


(152) 


Expanding in inverse powers of c^, we re-write the above equation of the 
stream-lines (I152|) as m- 


pA) 


dV{t) 

dt 


+ VW(x) 


Vp{x] 


I- 



V(t) 


dV(f) 

dt 


- 3V(t) • VIU ) V(t) + (21U + |V(t)p) VIU 


+ P + 3VW) + (4W - |V(t)|2]) Vp + (v{t) • Vp + V(t)| + O (i) 


(153) 
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The right hand side contains all the relativistic corrections to the Euler equations of 
Newtonian fluid dynamics m, d- 


7 Generalizations to complicated materials and curvilin¬ 
ear coordinates 


7.1 Perfect fluid plasma 

In this section we maintain field equations (ll()5al - [l()5ej) for the fields 4>(i)ap{x) 
4‘{o)apix)- We also retain equations (I88aj) and (I88cl) for the exterior metric 
and the normalized vector u^. Moreover, we stipulate the same equation as H1()7B for 
the interacting energy momentum stress tensor 0"^(x). 

A charged fluid or plasma satisfies the electromagnetic field equations (I5f)al 
I55c|) . The energy momentum stress tensor for this system is given by 


rpa/3 
^ (nvp) 


(x) 


p{x) + 



u°‘{x)u^{x) + p{x)g^QS^{x) + 0"^(x) + M°‘^{x), 


(154) 


where the subscript (nvp) denotes “non-viscous plasma”. The quantity A4“^(x) is 
furnished by (I58al) . Therefore, equation (|105dl) implies that 


„ d 


0 = 




= u 


dx!^ - 


0 + f 


+ (p+ ( pg?J^. ) -I- —Fg^. 


d 


/9 




dx0 V^^(o) 


jP 


(155) 

Here, we have used equation (11071) for 

Multiplying CSSl) by Ua (and summing), we derive the plasma continuity 
equation: 


A 

dxP - 




Ua 


d 


0/9 


dxP V^^(O) 


+ —F, 




(156) 


Substituting this last equation into (USSl), we deduce the generalized Euler equation 


(^’+ 3 ) +(n+^ 


d / 


dxP 


P 9 i^^ ]+-Fr 


= 0 . 


(157) 


From (ITTTl) . the equations for a stream line emerges as 


, P 

P + -2 
c"'J 


FX^{s) f a 


ds^ 


+ 


(3 7 


dXP{s) dX^{s) 


(0) 


ds ds 


'a 9{0)g.-j{x) dX>^{s) dX°^{s) 
(? ds ds 


\x°‘=X°^{s) 

(p{2>)<(^)) + 


I- 

(158) 
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One may add viscosity to the above system. For such a fluid, the energy 
momentum stress tensor is 




x) + 


u°‘{x)u^{x) 


du^{a 


dx° 


+ 






(9n"(x) 

dx^ 


+ 


-r]{x) - C(x) 


' af), . , U^{x)ul^{x) 
5 ( 0 ) (^) + - ^2 - 


du^ 


dx° 


(159) 


Here, r/(x) and (^{x) represent the shear viscosity and bulk viscosity respectively. Using 
we obtain, from (CMl) the following equation 


5r“5(x) _dT^£^){x) d 


0 = 


dxh 


dxh 


dxh 


5 ( 0 ) 


+ 


du^ 

dx^ 


+ ( 5(0) + 


j3fj \ (9n" 


dx° 


+ c-,. uro) + 




dx° 


(160) 


Here we have again stipulated the equation (IMI). 

Multiplying (Cnni) by Ua, we derive the continuity equation: 


d 

dxh L 


P + 




Un 


d 


afS 


dxh M^(o) 


jh 

+ -F,- 


d 

dxh 


5 ( 0 ) 


+ 


du^ 


+ 




5 ( 0 ) 


+ 


dx° 


+ (c - li) (sCol + 






dx^ 


(161) 


Substituting (dnu) into (Tmnii . we deduce the generalized Navier-Stokes equation: 

J5 


r 5 

( „„'y0\ 

dxh 

[P9io)) 


-F, 


+ 


d 

dxh 


4) + 


u'^u^\ du^ 


dx° 


+ ( 5(0) + 


0 

\ du'^ 


dx° 




+ 


u ' 


0u^\ du^ 


dx° 


( 162 ) 


7.2 Curvilinear coordinates and orthonormal frames 

Now we shall introduce curvilinear spacetime coordinates by transformation equations: 


x“ =X“(x), 
d (x^, x^, x^, x^) , 

—:—:—:— L ^ o 

d(x^, x^, x^, x"^) ’ 

x“ =X“(x). 


(163) 
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Here, we have assumed that the functions X°^ are of class The transformation 
rules for tensor fields are furnished by HH 


fa,.. dx-{x)dxf^{x) 

g^5 ■■ g^,, g^u P^.A )■ 1 ) 

The coordinate transformation (CMl) generate a new metric tensor field as 


follows: 


dapi^) ■ — d^i/ 


dX>^{x) dX’^(x) 
dxd 


dx°^ 
1-1 


:=[gai3{x)Y 

The Christoffel symbols are given by (compare with (ISIjl i 


P 7 } 


dg-faix) dgapix) dg,^{x) 


(165) 


dxd ' dx"^ dx° 

The covariant derivatives are defined by [TT| . pH] (compare with (|86al - l86Hl) l 


(166) 




If we replace the various tensor fields T'^^"pa..{x) by T°^^'/lu..ix) (in (I164jl i. and 
VqT’^'^'pct(x) (in (11671) . then all the constitutive equations are expressed in curvilinear 
spacetime coordinates 

However, for the sake of applications, we restrict ourselves to spatial curvi¬ 
linear coordinates only. In that case the following restriction is placed on equations 

(ITHHl) 

X* = X*(x), = X^(x, x^) := x^, X* = X*(x). (168) 

The tensor transformation rules (tTTlHl) then simplify in the obvious way: 


fb4.. ^4^ _ dXYx) dX^jx) dX'^jx) dX^{±) ^^4 

kli.. ’ g^a g^b g^k gf 

f^(x,x^) =T\{x,x^), 

and the metric tensor field in (tTHTHl reduces to 

dX^{±) dX\^) 


mnA.. 1^1 ^ )! 


(169) 


ft4(x,X^) =0, 

544(x,x‘^) =-1. 


9x* 9x1 ’ 


(170) 
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The Christoffel symbols in (|166l) boil down to 


^ dgqjjx) 9ffjfc(x) 

\ j k ) 2 dx^ dx^ dx°‘ 



(171) 


The laws of covariant differentiation (Unzi) imply that 




4 ^ _(9r“"^^4(x,x^) ^ f a 


dx^ 


+ 


j k 


b4 


ix,x^) - 


k 

j b 




4, _af“\4(i,x4) _ 9r“\4(x,x4) 


V 4 r “\ 4 ( x , X ^) = 


dx^ 


dx^ 


(172) 


If we now replace various tensor fields occurring in the 

constitutive equations (|154l - ITT)^ by , T “^"4 , T '^"^44 and (x, x^), 

-^T°'^"cd (X) 6tc., then we have converted all relevant equations into spatial curvi¬ 

linear coordinates. 

In physical applications, usually orthonormal or physical components of a 
tensor are necessary m- For that purpose we introduce three orthonormal vectors, 
A/i(x) A G {1,2,3}, in space. These vectors satisfy the orthonormality conditions: 


ffii(x)A}i(x)A^(x) = 6ab- (173) 

We define the inverse entries by 

[hf] ■■= [Ay-\ 

= 6 ), (174) 


(Here, the summation convention is also followed for capital roman indices.) 

By (|B73|) 3)11(1 (tmi we obtain 

Qijiic) = (lAShf (x)/x^(x). (175) 

The tensor transformation rules lead to 

^’^^CD..(x,x^) = Aia (x)l^f (x)..Ac(i)A}5(x)..f“'’-ij (x,x4). (176) 

Instead of Christoffel symbols, we require Ricci rotation coefficients |23] for 
the connexion. These are defined by: 


7abc(x) ■.=gji{x) (VfcA(4^ A]5 (x)Ac(x) = -'jBAci^), 

7"bc(x) :=5"^‘^7JBc(x) = 'jABcix). (177) 
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The appropriate covariant derivatives can be characterized by: 


VjT 


'A4 


B4(x,x^) := A*j(x)^r 


iA4 


(i, - 7^^j(i)r%4 + 7^5j(i)r 


.n 


B4 


nA4 


D4- 


(178) 


If we now replace the tensor helds appearing in the 

constitutive equations (llbdl - ITT)^ by (x, x^), T "^^^4 (x,x^), V (x,x^), 

-^j'AB^^ (x, x^) etc., then we have transformed all the required equations into the 
orthonormal or physical frame. Physical measurements correspond to quantities ex¬ 
pressed in this frame. 

Before concluding, we shall now explore a special example which is most useful 
for continuum mechanics. This example involves orthogonal coordinate systems in 
Euclidean three-space. The equations jnni) reduce to 


[5ii(x)] 


[/ii(x)]^ 0 0 

0 [/i2(x)]^ 0 , 

0 0 [/i3(x)]^ . 


(179) 


/ii(x) > 0, Y^det [gij] = /ii(x)/i 2 (x)/i 3 (x) > 0. 

The non-zero Christoffel symbols, from (dZIl) and (dzni) are summarized by: 


1 

1 1 
1 

1 2 


d 

dx^ 

d 

dx"^ 


In hi, etc. ; 
In hi, etc. ; 


(180) 


1 

2 2 


— 2 ^ 

=-h 2 [hi] ^^lnh 2 , etc. 


The orthonormal (or physical) vector components from (I178B and (11791) are 
given by 

A*^(x) = [hi{±)]~^6\, /x^(x) = hi{±)6j , (181) 

and the non-zero Ricci rotation coefficients from (EZZI), (CHOI) and dHU) are furnished 
by: 


7ABC ■=^AE^%c^ 
7(2)(1)(2) 

7(1)(2)(1) = “ 
7(1)(3)(1) = “ 


1(A){B){C) '■= lABC = 1BAC-, 



d 

In hs 

2, 7(3)(1)(3) = 

“9x1 


d 

In hs 

1, 7(3)(2)(3) = 

9x2 


9 

In h 2 

1, 7(2)(3)(2) = 

9x3 


(182) 
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Also, the gradient of a scalar field is given by 


= {Ha) ^ 

V(i)<^(i) =(hi)“^^0(x) # Vi(^(i). 

In (11 there is no summation over A. 

The divergence of a vector field is furnished by 


ViT{iL) = VaT^{^) = {hih2h3)-^ \ ^ 


(hlh2/l3)r*(x) 


(183) 


(184) 


The curl is given by: 


V X A 


ijk 


^ j-^k ^ k 


1 


ijk 


V X A 


1 B 


2 v^det[5afe] 1 " ■'J 2(/ii/i2/i3) 

(io(x)) - A- « (ic(x)) 


dAk{±) 

dyj 


2^ 


BCD 


‘f 


BCD 


V X A 


11 


V X A 


{he) 

1(1) 


-1 


dx^ 


[M^))-{hD) i^(ic(i)) 


dAj (x) 
3x^ ’ 


(185) 


Finally, the Laplacian operator is furnished by 


V2iT(x) ;= 5*iViVjTT(x) = (5 ^'®VaVbIT(x) 


1 

a 

[ hih^ 

aiy(x) 

w ^ 

( h^hi 

aiy(x)\ 

(/11/12/13) 

a^i 

hi 

dx^ 

1 dx^ 

^ h 2 

dx"^ J 


(186) 

d f hi h 2 5IT(x)\ 

dx^ y hs dx^ j 
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